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CAUSALITY  DIFFICULTIES  WHICH  CAN  ARISE  IN  MODELLING  THE  REFLECTION 
OF  A  NORMALLY  INCIDENT  WAVE  FROM  A  LOSSY  PLANAR  SURFACE 


INTRODUCTION 

The  reflection  of  an  acoustic  wave  normally  incident  on  a  lossy  planar 
surface  is  of  theoretical  interest  owing  to  its  application  to  echo-reduction 
problems.  Frequently,  the  loss  mechanism  of  the  material  of  the  surface  is 
treated  theoretically  by  generalizing  the  solution  for  the  lossless  case  via 
introduction  of  complex  wavenumbers,  although  this  treatment  is  only 
approximately  correct  [1].  If  this  approach  is  used,  and  insufficient  care  is 
exercised  in  the  choice  of  the  frequency  dependence  of  the  material  properties 
of  the  substance  comprising  the  lossy  surface,  then  the  causality  condition 
may  not  be  satisfied. 

This  possibility  is  investigated  in  the  present  work  by  analytically 
computing  the  reflected  pressure  for  three  simple  material  models.  The 
material  of  the  planar  surface  is  modelled  to  have  a  frequency-independent 
phase  speed  in  all  three  cases.  The  lossy  nature  of  the  surface  is  chosen  to 
depend  on  frequency  such  that  it  (1)  is  constant,  (2)  varies  linearly  with 
frequency,  or  (3)  varies  quadratically  with  frequency.  The  approach  to  the 
reflection  problem  is  to  decompose  the  incident  pressure  field  into  its 
Fourier  plane-wave  components.  The  interaction  of  each  of  these  components 
with  the  surface  is  then  evaluated  using  the  well-known  result  for  the 
reflection  of  a  plane  wave  from  a  semi-infinite  material  region  (allowing  the 
wave  velocity  in  the  material  and  Its  associated  wave  vector  to  be  complex). 

It  will  be  shown  that  when  this  approach  is  applied  to  the  three  material 
models  described  above,  only  case  (1)  (frequency-independent  loss)  results  in 
a  reflected  pressure  with  appropriate  causality  properties.  Hence,  caution 
must  be  exercised  if  model  (2)  or  (3)  is  applied  to  situations  in  which  the 
attendant  causality  difficulties  would  be  unacceptable.  (Such  causality 
difficulties  can  arise  when  a  truncated  Fourier  solution  is  used,  due  to 
Gibb’s  phenomenon.  This  is  not  the  case  In  what  follows,  however,  since  the 
required  integrals  are  evaluated  exactly  by  analytical  means.) 

This  problem  has  been  considered  previously  [2];  however,  approximations 
were  used  to  evaluate  the  relevant  integrals.  These  approximations  were  not 
necessary  in  the  present  research.  Also,  although  it  is  clear  from  Fig.  3  of 
Ref.  2  that  causality  difficulties  arise  from  these  models,  no  specific 
discussion  of  this  is  given  in  that  reference.  The  purpose  of  the  present 
work  is  to  identify  these  causality  difficulties,  explain  why  they  arise,  and 
discuss  their  resolution. 


PRELIMINARY  CONSIDERATIONS 

The  normally  incident  pressure  field  is  assumed  to  be  of  the  form 
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pi(t)  =  Vo 
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The  exact  solution  for  the  reflected  field  is  quite  complicated  to  obtain, 
involving  boundary  conditions  that  not  only  specify  the  usual  displacements 
and  stresses  but  also  involve  the  acoustic  temperature  [3].  However,  there  is 
a  large  body  of  empirical  evidence  that  demonstrates  that  acoustic  losses  in 
materials  can  be  modelled  fairly  accurately  in  many  cases  by  simply 
generalizing  the  solution  for  the  lossless  case  by  suitably  allowing  the 
wavenumber  (and  the  wave  speed)  to  be  complex.  This  latter  approach  is 
adopted  here.  Thus,  each  of  the  plane-wave  components  of  the  pressure  field 
within  the  lossy  material  itself  is  assumed  to  be  representable  by  functions 
of  the  form 


pt(w) 
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where  k  =  00/(0^)  ~  complex  wavenumber 


frequency,  c 


ph 


Here ,  oo  is  the  angular 
is  the  phase  speed  in  the  material,  and  a  is  the  loss  factor. 


The  Fourier  components  of  the  incident  field  are  computed  from  p^(t)  using 


Ai(u)) 
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and,  due  to  the  reality  of  the  incident  field,  the  inverse  Fourier  transform 
may  be  written 


Pt(t)  =  2 Re 


00 

f  A^((jo)e'*'a)tda)  . 
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(4) 


The  symbol  r(w)  is  used  to  represent  the  (complex)  reflection  coefficient, 
where  the  well-known  result  [4],  generalized  using  a  complex  velocity,  is 


r(m) 
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Here,  pQ  and  cQ  are  the  density  and  phase  speed  of  the  fluid  (assumed 
lossless)  in  contact  with  the  lossy  medium.  Also,  and  Cj-  are  the  density 
and  complex  sound  velocity  in  the  lossy  medium  itself,  with  =*  0)/k  .  The 
reflected  pressure  is  computed  from  r(u>)  using  the  formula 

oo 

Pr(t)  =  2 Re  J  r(0))Ai(a))ela)tda)  .  (6) 
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ANALYTICAL  COMPUTATION  OF  THE  REFLECTED  PRESSURE 

Before  proceeding  further,  it  is  helpful  to  note  that  r(co)  can  be  put 
into  a  convenient  form  with  the  aid  of  the  following  definition: 


,  PoCo 

a+  ■  pi  i  — 


•ph 


Using  this,  r(u))  becomes 


r(to) 


a  (0  +  ia  p  c 
- _ o  o 

a .  a)  -  ia  p  c 
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(7) 


(8) 


The  reflected  pressure  wilL  next  be  evaluated  analytically  for  each  of  the 
three  simple  models  of  interest. 


Constant  Loss  Factor 

When  Eq.  (8)  is  substituted  into  Eq.  (6)  (with  a  =  constant),  the 
resulting  integral  can  be  put  into  a  form  amenable  to  the.  use  of  standard 
integral  tables  [5]  via  partial  fraction  analysis.  However,  the  relevent 
tabular  expressions  involve  functions  possessing  a  branch  cut,  which  can  be 
particularly  troublesome  if  caution  is  not  exercised.  Following  is  a  list  of 
integrals  that  result  from  a  careful  analysis  of  this  branch  cut: 
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ci(ag)sin(a8) 

-cos(a8)si(a8) 
where  Re(8)>o  and  Im  (8)>0, 
but  not  simultaneously  =  0. 
Also,  a  is  real  and  >  0. 
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cos(ax)dx 
x  +  B 


-sin(a8)si(a8)-cos(a8)ci(a0) 
where  Re8>o  and  Im8>0, 
but  not  simultaneously  =  0. 
Again,  a  is  real  and  >  0. 


(10) 
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sin(ax)dx 
Y  -  x 


[ci(ay)  -  ir]sin(ay) 
-cos(ay) [si(ay)  +  "]  , 
where  Re(Y)>0,  and 
Im  y  *  0.  (If  Im  y  =  0, 
itr  term  is  dropped.) 
Again,  a  is  real  and  >0. 
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+  cos(ay) [ci(ay)  -  irr ]  , 
where  Re(y)>0,  and 
Im  yW  (If  Im  y  =  0,  iir 
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a  is  real  and  >0. 


(12) 


In  Eqs .  (9)  thru  (12),  ci  and  si  are  the  cosine  integral  and  sine  integral, 
respectively  [6].  Using  these  results,  analysis  of  the  required  integral  is 
straightforward ,  though  tedious.  The  result  is: 


p  (t)  =  p  Re 
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It  should  be  noted  that  in  the  limit  as  a-K),  Eq.  (13)  becomes  the  well- 
known  lossless  result  pr(t)  =  (a_/a+)  p0sin(wot)  (1-e  ). 


Loss  Factor  Proportional  to  Frequency 

Next,  the  function  a  is  allowed  to  assume  the  form  a  =  a  (w/w  )  .  It 

should  be  noted  from  Eq.  (8)  that,  in  this  case,  the  reflection  coefficient 
r(w)  is  frequency-independent.  Hence,  this  function  may  simply  be  factored 
out  from  the  integral  of  Eq.  (6).  [In  what  follows,  the  symbol  rQ  is  used  to 
represent  the  complex  constant  r(w).]  However,  the  resulting  integral  does 
not  simply  reproduce  the  incident  pressure  p^(t)  [compare  Eq.  (4)].  This  is 

because  the  imaginary  part  of  the  integral,  which  does  not  contribute  to  the 
computation  of  pi(t),  does  contribute  to  the  reflected  pressure  due  to  the 
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effect  of  the  complex  (constant)  reflection  coefficient.  It  is  also 
interesting  to  note  that  since  r(o>)  -  rQ  is  constant  in  this  case,  it  follows 

from  Eq.  (5)  that  the  complex  sound  speed  c^  is  also  a  constant. 

The  integrals  arising  in  this  case  can  still  be  evaluated  using 
Eqs.  (9)  thru  (12).  Rather  than  directly  displaying  the  reflected  pressure  in 
this  case,  it  is  somewhat  more  convenient  to  introduce  a  new  quantity  I,  where 


oo 


(14) 


0 


It  should  be  noted  that  this  is  also  the  integral  appearing  in  Eq.  (4).  The 
reflected  pressure  then  becomes: 


Pr(t)  =  Re(rQI) 


(15) 


T  >  t  >  0  , 


where 


irn  t 
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e  iUot  [ci(u)Q  -  iY)t  +  isi(coQ  -  iY)t] 


i  T  '  ' 
+  —  Im ,  e 


TT 


(16) 
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The  real  part  of  Eq.  (16)  is  equal  to  pQ(l-e  )sin((Dot),  as  is  expected  based 


on  its  relationship  to  Eq.  (4),  the  Fourier  representation  of  the  incident 
field. 
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Loss  Factor  Proportional  to  Square  of  Frequency 

2 

Finally,  the  function  a  is  allowed  to  assume  the  form  a  »  a^m/w^)  . 

When  this  is  substituted  into  Eq.  (8),  we  obtain  an  expression  for  r ( oo)  very 
similar  to  that  obtained  for  the  case  a  =  constant  (in  the  present  case,  a 
factor  of  a)  is  present  in  the  imaginary  parts  of  the  numerator  and  the 
denominator  of  r(w),  whereas  in  the  case  a  =  constant,  a  factor  of  a)  appeared 
in  the  real  parts).  Hence,  a  very  similar  partial  fraction  decomposition  is 
possible,  and  Eqs.  (9)  thru  (12)  again  suffice  to  evaluate  the  resultant 
integrals.  Before  displaying  the  final  expression  for  the  reflected  pressure, 
a  few  definitions  are  helpful: 
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In  terms  of  these  quantities,  the  reflected  pressure  becomes: 
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P^(t)  -  p. 
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T  >  t  >  0  .  (17) 


CAUSALITY 


Consider  next  the  consequences  of  the  above  results  [Eqs.  (13),  (15)  and 
(17)]  in  regard  to  causality.  Although  no  difficulties  arise  for  Eq.  (13) 

(the  case  where  a  =  constant),  this  is  not  so  for  Eqs.  (15)  and  (17)  (where  a 
is  proportional  to  frequency  and  proportional  to  the  square  of  the  frequency 
respectively).  The  most  obvious  way  to  see  that  there  are  difficulties  in  the 
latter  cases  is  to  note  that  the  results  are  dependent  on  the  pulse  duration 
t,  a  though  the  results  presented  are  valid  for  times  less  than  t.  This  is 
clearly  inconsistent  with  causality  requirements,  since  it  means  that  the 
end  Ox  the  incident  pulse  influences  the  reflected  pressure  prior  to  its 
arrive1  at  the  reflecting  surface.  Although  not  as  obvious,  careful  analysis 
or  these  expressions  also  shows  that  when  aQ  £  0  then  pr(t)  #  0 


at  t  0.  Hence,  once  again,  these  expressions  contradict  the  causality 

xequirement  that  no  reflected  wave  should  arise  prior  to  the  arrival  of  the 
incident  wave. 


There  is  a  simple  alternative  way  to  see  that  causality  difficulties  will 
arise  for  these  models.  This  is  by  computing  the  impulse  response  of  the 
given  system  [i.e.,  by  computing  the  reflected  pressure  that  would  arise  for 

an  incident  field  of  the  form  ?i(t)  =  p0<5(t),  where  S(t)  is  the  Dirac  delta 

function].  This  is  done  by  computing  the  inverse  Fourier  transform  of  the 
function  r(u>).  As  an  example,  consider  the  model  in  which  a  is  linearly 
dependent  on  frequency.  Recall  that,  in  this  case,  r(u>)  is  a  complex  constant 
(chosen  to  be  rQ).  Thus 


r(t)  5  F  1  [r(rn)]  =  2Re 


/' 


loot 


dco 


(18) 


or 


r(t)  =  2Re{ro[7T<5(t)  +  i]}. 


(19) 


Hare,  F  ^  refers  to  the  inverse  Fourier  transform.  From  Eq.  (19) 
that  when  Im(rQ)  *  0,  then  r(t)  *  0  for  t  <  0.  Hence,  there  is  a 
response  prior  to  t  =  0  to  an  impulse  at  t  =  0. 


it  is  clear 
non-zero 
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Although  the  above  conclusions  may  be  disturbing,  they  are  not 
unreasonable  for  the  following  reason:  In  order  for  a  signal  to  be  physically 
realizable  (i.e.,  causal),  it  must  conform  to  specific  mathematical 
conditions.  In  fact,  the  real  and  imaginary  parts  of  the  frequency  spectrum 
of  the  signal  must  be  Hilbert  transforms  of  each  other  [7,8].  In  view  of  this 
rather  stringent  mathematical  requirement,  it  is  unlikely  that  the  reflected 


pressure  pr(t)  arising  from  the  interaction  of  a  wave  with  a  material  of 

arbitrarily  selected  material  properties  will  be  consistent  with  the  required 
transforms.  In  fact,  in  retrospect,  it  is  remarkable  that  the  case  a  = 
constant  does  not  result  in  causality  difficulties.  The  most  probable 
rectification  of  the  problems  arising  in  the  latter  two  cases  is  to  select  an 
appropriately  frequency-dependent  material  phase  speed  (i.e.  ,  it  must  be  that 


the  only  functional  form  for  a  which  is  consistent  with  cph  =  constant  is  a  = 
constant  as  well).  It  is  not  intended  here  to  imply  that  such  approximations 


as  c 


ph 


constant 


and  a  proportional  to  u> 


2 


are  not  reasonable  for  many 


materials.  It  is  only  intended  to  point  out  that  one  consequence  of  using 
these  approximate  expressions  is  the  breakdown  of  the  causality  condition. 
(For  guidance  in  the  proper  selection  of  the  frequency  dependence  of  the  real 
and  imaginary  parts  of  material  properties,  see  Refs.  1  and  8.) 
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